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ABSTRACT 

The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spec- 
troscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to ob- 
tain spectra for up to ~ 300 000 galaxies over 280 square degrees, to a limiting magnitude of 
r pct < 19.8 mag. The target galaxies are distributed over < z < 0.5 with a median redshift 
of z » 0.2, although the redshift distribution includes a small number of systems, primar- 
ily quasars, at higher redshifts, up to and beyond 2=1. The redshift accuracy ranges from 
a v w 50kms _1 to a v k, 100 kms -1 depending on the signal-to-noise of the spectrum. Here 
we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps 
involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated 
one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range 
of 3750 < A < 8850 A at a resolution of R m 1300. The final flux calibration is typically 
accurate to 10 — 20%, although the reliability is worse at the extreme wavelength ends, and 
poorer in the blue than the red. We present details of the measurement of emission and absorp- 
tion features in the GAMA spectra. These measurements are characterised through a variety 
of quality control analyses detailing the robustness and reliability of the measurements. We 
illustrate the quality of the measurements with a brief exploration of elementary emission line 
properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence 
of the Balmer decrement, consistent with previously published results, and explore further 
how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass 
and redshift dependencies of the [NII]/Ha vs [OIII]/H/3 spectral diagnostic diagram, com- 
monly used to discriminate between star forming and nuclear activity in galaxies. 

Key words: galaxies: evolution - galaxies: formation - galaxies: general 



1 INTRODUCTION 

Galaxy surveys that have moderate resolution (R ~ 1000 — 2000) 
optical spectroscopic measure ments, such as the T wo-degree Field 
Galaxy Redshift Survey (2 dF: IColless et al.ll200l1) and the Sloan 



separa tion, and target prioritisation is presented by iBaldrv et al 
(2010), with the tiling algorithm described by Robotha m et al 
( |2010|) and the photometric analysis by Iffill et alj d201ll) . Stel- 



Digital Sky Survey (SDSS; York et al. 2000) are among the most 



productive resources for understanding galaxy formation and evo- 
lution. The astrophysical information encoded in the rest-frame op- 
tical region of the spectrum is among the most well-understood 
and well-calibrated aspect of galaxy evolution studies, and pro- 
vides a wealth of detail regarding the physical processes occurring 
within galaxies. This ranges from quantitative measurements of star 
formation rate (SFR), metallicity, velocity dispersion, obscuration, 
and more, through to diagnostics distinguishing between the pres- 
ence of star formation or an accreting central supermassive black 
hole (an active galactic nucleus, AGN). In combination with broad- 
band photometric measurements spanning ultraviolet through to ra- 
dio wavelengths, the redshift and other physical information from 
galaxy spectra provides a powerful tool for exploring the details of 
galaxy evolution. 

The Galaxy And Mass Assembly (GAMaQ survey is a 
large multiwavelength photo metric and spectroscopic galaxy sur- 
vey ( iDriveret all 120091, |20 111) that provides exactly this compre- 
hensive selection of photometric and spectroscopic data. The key 
scientific goals are to use the galaxy distribution to conduct a 
series of tests of the cold dark matter (CDM) paradigm, in ad- 
dition to carrying out detailed studies of the internal structure 
and evolution of the galaxies themselves. The scientific moti- 
vation for GAMA, the survey footprint, data proc essing, cata- 
logue construction and quality control are described bv lDriver et al.l 
d201ll) . The target selection, including survey masks, star-galaxy 
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lar masses for th e GAMA galaxies have been quantified by 
iTavlor et al.l d201 1[). and the lo w redshift stellar mass function is 
detailed in Bald rv et alj d2012h. The broadb and luminosity func- 
tions are derived by lLovedav et al. I d2012l) . and the Ha lumi- 
nosity functions and evolution presented in iGunawardhana etHEl 
d2013l). Gal axy nebular metallicity measuremen ts are detailed in 
iFoster et"al] d2012l) and lLara-Lopez et al.l d2013l). Ga l axy groups 
in GAMA have been quantified by iRobotham et"al] 1201 ll) . and 
galaxy structural parameters measured by Kelvin et al. Id2012j). 



The GAMA survey used 68 nights of observing time on the 
Anglo- Australian Telescope (AAT) over 2008-2010. This time was 
used to conduct a highly complete survey in three Equatorial fields 
to r pc t < 19.4 mag over a total of 144 deg 2 , 48 deg 2 of which was 
observed to the deeper limit of r pct < 19.8 mag. This initial phase 
of the survey, usually referred to as "GAMA I," allowed the acqui- 
sition of over 112 000 new galaxy spectra and redshifts, for a total 
of over 130 000 redshifts in the original GAMA survey area. Subse- 
quently the survey has been extended, with the award of 1 10 nights 
of AAT time over 2010-2012, to expand the survey by including 
two Southern fields and broadening the three Equatorial fields (re- 
ferred to as "GAMA II"). This expands the total survey area to 
280 deg 2 , while achieving a uniform depth of r pct < 19.8 mag over 
the full survey region. The goal is to compile ~ 300 000 galaxy 
spectra over this area. At the time of writing, we have already ob- 
tained over 220 000 spectra. In detail, to date GAMA has observed 
224465 spectra of which 222 294 are galaxy targets. These num- 
bers include repeat observations, and not all are main survey tar- 
gets, as they include "filler" targets that take advantage of fibres 
unable to be allo cated to main surv ey targets on any given obser- 
vation plate (see Bald rv et alj|2010l) . Including spectra from other 
surveys within the GAMA regions, (SDSS, 2dFGRS and others, see 
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Baldry et ai] |20ld. for details), we have 299 980 spectra, of which 
297 067 are of galaxy targets (not all are main survey targets). We 
have 233 777 unique galaxy tar gets, of which 215 458 (92.2%) have 
redshift quality nQ ^ 3 (see lDriveretaill201ll for definition of 
nQ, but in brief nQ — 3 or nQ = 4 correspond to reliable red- 
shifts). 

The survey has already led to a number of published re- 
sults making use of the detailed emission and absorption line 
measurements from the GAMA spectroscopic data. These include 
identification of the lowest-mass star forming galaxy population 
l lBroughetai]|201lh . a self-consistent approach t o galaxy star for- 
mation rate estimates and the role of obscuration (Wiiesinghe et al. 
12011 al ibi), and evidence for a star formation rate dependence 
in the high-mass slope of the stellar initial mass function 
dGunaw ardhana et al.l201ll) . among many other GAMA team pub- 
lications (for a full and current list see the team web page), along 
with additional work led by collaborating surveys such as Herschel- 
ATLA50. 

In order to best facilitate subsequent scientific analyses of pub- 
lic survey data, it is crucial to provide full details of the observa- 
tions, processing and data product derivations (e.g., iBolton et al.l 
|2012|) . Here we describe the spectroscopic pipeline process- 
ing for the GAMA survey. This encompasses an overview of 
the observations (§[2}, the steps involved in processing the raw 
two-dimensional spectroscopic images through to extracted one- 
dimensional spectra and the initial redshift measurement process 
(§[3}, and flux calibration of the one-dimensional spectra (§|4). We 
also present the processes used in measuring the emission and ab- 
sorption features of the spectra (§|5]( that are recorded in the GAMA 
database and made publicly available through the staged data re- 
leases. Note that the GAMA spectroscopic reduction and analysis 
pipeline is still evolving as we continue to improve some aspects. 
Here we describe the pipeline that was used to construct the final 
GAMA I dataset. It was this dataset that has been used in the in the 
various investigations cited above. The spectra and associated mea- 
surements that will be available in GAMA DR2, the public data 
release due in January 2013, also rely on the pipeline described 
here. For the purposes of this paper we use the data associated with 
GAMA SpecCat v08. 

Throughout, all magnitudes are given in the AB system, and 
we assume a cosmology with Ho =70kms~ 1 Mpc~\ Q.m = 0.3 
and Q,a = 0.7. 



2 OBSERVATIONS AT THE AAT 

The G AMA spectroscopic observations use the AAOmega spectro- 
graph dSaunders et aill2004 ISmith et al.ll2004l Isharp et al.l 120061) 
on the 3.9 m Anglo-Australian Telescope (Siding Spring Obser- 
vatory, NSW, Australia) for measuring the spectra of the target 
galaxies. This spectrograph is stationed in the thermally stable en- 
vironment of one of the telescope Coude rooms. AAOmega pos- 
sesses a dual beam system which allows coverage of the wave- 
length range from 3750 A to 8850 A with the 5700 A dichroic used 
by GAMA, in a single observation. Each arm of the AAOmega 
system is equipped with a 2k x 4k E2V CCD detector and an AA02 
CCD controller. The blue arm CCD is thinned for improved blue 
response. The red arm CCD is a low fringing type. The grating 
used in the blue arm is the 580V, centred at 4800 A, which has a 

2 |http : / /www. h-atlas ■ org/| 



dispersion of 1 A/pixel and gives a coverage of 2100 A. The 385R 
grating is used in the red arm, centred at 7250 A. This grating has 
a dispersion of 1.6 A/pixel and gives a coverage of 3200 A. This 
leads to spectra with a resolution that varies as a function of wave- 
length, from R fa 1000 at the blue end up to R w 1600 at the red 

end. 

The "Two-degree Field" (2dF) instrument dLewis et al 
consists of a wide field corrector, an atmospheric dispersion com- 
pensator (ADC), and a robot gantry which positions optical fibres 
to an accuracy of 0"3 on the sky. T he fibres have a 2" diameter pro- 
jected on the sky dLewis et al.l(2002l) . A tumbling mechanism with 
two field plates allows the next field to be configured while the 
current field is being observed. The 392 target fibres from 2dF are 
fed to the AAOmega spectrograph, and eight guide fibre-bundles 
are used to ensure accurate telescope positioning over the course of 
each exposure. For GAMA observations, the 2dF robot fibre posi- 
tioner is used to configure typically 345 fibres to observe galaxies 
within a two degree field on the sky. Due to a varying number of 
damaged or unusable fibres (typically around 20), the actual num- 
ber of fibres able to be used for galaxy targets is not constant. To 
quantify this, the first-quartile/median/third-quartile number of fi- 
bres on galaxy targets for GAMA I was 332/345/348. For the full 
survey to date, these numbers are 324/341/348. This distribution 
has remained fairly steady over the duration of the survey so far. 
Around 25 additional fibres are used to measure the sky spectrum in 
ea ch field. Sky positio ns were identified using a sky mask, detailed 
in iBaldrv et alj d20ld) . Three fibres are allocated to spectroscopic 
standard stars. 

The integration time for each GAMA field is typically 60 min, 
split into three 1200 s exposures. Accounting for the read-out time 
of the CCDs (2 min) and the acquisition of the calibration frames 
comprising flat-fields and arc-line exposures for wavelength cali- 
bration, the time spent on each field is well-matched to the time re- 
quired for the 2dF positioner to configure the following observing 
plate, ensuring an efficient overall survey strategy. Between 20 and 
30 bias frames are taken during each observing session. Starting in 
2011, we also began to use dark frames to refine the calibration, 
with from 10 to 30 dark frames, each of 1200 s exposure, being 
acquired each observing session. 



3 DATA PROCESSING AND REDSHIFT 
MEASUREMENTS 

3.1 Obtaining ID spectra 

The raw data ar e processed using software developed at the AAO 
called 2DFDR, dCroom et al] 12004 Isharp & Birchalll |2010|) . The 
2DFDR processing applies the standard sequence of tasks for ID 
spectral extraction from 2D images. This includes bias subtraction, 
flat-fielding, fibre trace (or "tramline") fitting, and wavelength cal- 
ibration. First, a master bias (and for GAMA II data a master dark) 
are created using the available bias (and dark) frames. For each 
new plate configuration observed, the raw AAOmega frames are 
run through 2DFDR at the telescope to provide the processed spec- 
tra. 

The standard parameters are used in running 2DFDR, with the 
following modifications. We consider not only the master bias (and 
master dark), but also an overscan correction using a ninth- order 
polynomial fit for the blue spectra, and second-order for the red 
spectra. The high order is necessary in the blue due to a very strong 
gradient in the first ~ 100 pixels that is not well-modelled by lower 
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Figure 1. Example GAMA spectra for a selection of galaxies, illustrating high quality spectra and spectra affected by instrumental and processing artifacts. The 
spectra have been smoothed with a five-pixel running boxcar average to aid in clarity of display. From top to bottom: Star forming galaxy spectrum; Absorption 
line galaxy spectrum; Spectrum affected by fringing; Spectrum affected by bad splicing. Each spectrum includes an inset showing the S PSS colour galaxy 
image, as well as identifying common emission or absorption features. Each galaxy has its GAMA ID, redshift and stellar mass (from iTavlor et alfeoi II) 
listed, along with Balmer decrement, Ha luminosity and star formation rate for star forming galaxy spectra. Only about 3% of GAMA spectra are affected by 
fringing or bad splicing. 
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order fits. The wavelength solution is determined by a third- order 
polynomial fit to the arc-lines, with the solution tweaked using a 
first-order (blue) or third-order (red) polynomial to the sky lines. 
The higher order is required for the red CCD as there are more 
sky lines and the sky is brighter at the red end of the spectrum. 
The wavelength calibration is referenced to arc line wavelengths 
in air rather than vacuum, contrary to the convention adopted by 
SDSS, and is accurate to better than 0.1 A, as measured from key 
strong sky line features. Cosmic rays are cleaned in each object 
fram e using an implem entation of Laplacian cosmic ray identifica- 
tion ( Ivan Do kkum 2001), and applying clipping thresholds of 10 a 
in the blue and 5 a in the red. The scattered light is subtracted 
assuming a first-order polynomial fit. The throughput calibration 
method considers a flux weighted value of the night sky emission 
lines to normalise the fibre throughputs. 

Extraction of spectra is performed by first identifying fibres 
in a flat field frame and then fitting their locations as a function 
of CCD position using a model of the spectrograph optical distor- 
tion. Once the tramlin es are located, a minimum variance G aussian 
weighted extraction JSharp & Birchaiill2010l; |Hornel fl986) is used 
to obtain the flux in each fibre per spectral pixel. While optimally 
weighted, this does not take into account cross-talk between fibres, 
but given the restricted dynamic range of the GAMA targets (a 
range of less than three magnitudes) the level of cross-talk has been 
shown to be negligible ( Sharp & Birchalll2010L their Figure 2). For 
each ID spectrum there exists a variance array determined from 
the Poisson noise in the bias corrected 2D frame, the read noise 
and the gain, and propagated through the reduction pipeline. Ex- 
amination of repeat spectra shows that the uncertainties per pixel 
are well characterised by the measurements in the variance array, 
except around the strong 5577 A sky line, and at the extreme ends 
of the wavelength range. Where there are differences between the 
repeat measurements and the variance arrays due to such systemat- 
ics, the differences are always less than a factor of 1.4. 

The initial sky subtraction is performed using the ~25 fibres 
allocated in each plate to sky positions. A combined sky spectrum 
is made by taking the median of the corresponding pixels in each 
of the normalised sky fibres, discarding the two brightest sky fibres 
to avoid potential problems in the event of inadvertent non-sky flux 
(due to an asteroid, passing satellite, or other moving object per- 
haps), at the sky fibre location. The continuum sky subtraction ac- 
curacy is typically 2-3% of the sky level, although for especially 
strong sky lines such as that at 5577 A, th e residuals can be worse 
(for details see Sh arp & Parkinson! |2010|) . This is then followed 
by an improvement to the sky subtraction base d on subtracting a 
comb ination of principal component templates l lSharp & Parkinson] 
|2010|) . This reduces the amplitude of the sky subtraction emission 
line residuals to below 1% in most cases. 

Strong atmospheric telluric absorption features in the red part 
of the spectrum need to be corrected for. The telluric correction in- 
volves constructing a flux and variance weighted combination of 
all the spectra in a given field, which is then iteratively clipped to 
remove residual emission or absorption features (such as galaxy 
emission lines). This process relies on the fact that in any one 
field there are a broad range of galaxy redshifts so features are not 
present at the same wavelength in many spectra. The resulting aver- 
age spectrum is fit by a low order polynomial in the region around 
the telluric features, while excluding the regions where the absorp- 
tion is present. Dividing through by this polynomial fit results in 
a telluric correction spectrum which is set to be equal to unity ev- 
erywhere outside of the telluric absorption bands. The correction 



is then implemented by dividing each individual spectrum by this 
correction spectrum. 

Finally, 2DFDR splices together the blue and red spectra for 
each galaxy by doing a first-pass flux calibration to best match the 
spectra at the splice-wavelength (5700 A). The pixel scale is 1 A in 
the blue and 1.6 A in the red, although during the splicing step the 
red spectra are resampled to the same pixel scale as the blue, so 
the final pixel scale is ~ 1 Apix -1 . This is done with a quadratic 
interpolation ensuring that flux is conserved, and with appropriate 
treatment of masked or otherwise bad values. The same resampling 
is applied to the variance arrays to correctly propagate the errors. 
The overlap region between the blue and red spectra is 250 A, from 
typically 5650 A to 5900 A, although this varies slightly from spec- 
trum to spectrum depending on the location on the detector. 

As 2DFDR is continuing to be developed and improved, and to 
mitigate against reduction mistakes at the telescope during observ- 
ing, the entire GAMA dataset is periodically re-reduced, to ensure 
that the final spectroscopic data products are homogeneous and in- 
ternally self-consistent. The version of 2DFDR used in producing 
the final GAMA I spectra was 2DFDR v4.42. 

Fig. shows examples of GAMA spectra after the ID ex- 
traction and flux-calibration process (see §|4]below) are completed. 
This Figure shows two high quality spectra, and two poor quality 
spectra illustrating some of the instrumental and processing limi- 
tations in the survey. The two high quality spectra are an emission 
line object and an ab sorption line obje ct, both with redshift qual- 
ity of nQ — 4 (see Driver et al.ll201 it for the definition of red- 
shift quality flags and conventions). The first poor quality spectrum 
shown gives an example of fringing, visible as the high-frequency 
oscillation in the continuu m level, and acc ompanied by poor re- 
moval of the sky features JSharDetal.ll2006ll2013h . The fringing, 
which is time-variant, is only present for some fibres, and arises 
due to air gaps in the glue between the prism and the ferrule. Over 
time these fibres have been re-terminated with new glue and new 
ferrules. While it doesn't resolve the problem for existing GAMA 
spectra that are affected, the AAO has recently completed a total re- 
placement of all 2dF fibres with optimal glue and ferrules that has 
now eliminated this problem. In this particular example spectrum 
a redshift is still able to be reliably measured, with nQ = 4. The 
second poor quality spectrum shown is an example of a bad splice, 
characterised by a dramatic change in continuum level at the splice 
wavelength of 5700 A. This feature is a consequence of poor con- 
tinuum level estimation due to poor flat-fielding in one or both of 
the blue and red arms, in estimating how to scale the two compo- 
nents for the splice. The reliability of the splicing of AAOmega 
spectra is an area that is the subject of ongoing work, both at the 
AAO with continued development of 2DFDR, and within GAMA 
through investigation of independent flux calibration processes for 
the blue and red arms separately. For this spectrum, again a redshift 
is still able to be measured with high reliability, nQ — 4. These 
examples illustrate a key point, that a high quality redshift mea- 
surement can be obtained from a poor quality spectrum, although 
subsequent measurements of emission or absorption features for 
that spectrum may not be reliable. 

3.2 Redshift measurements 

Redshifts are measured from the one-dimensional galaxy spectra 
as soon as each field is fully reduced using the above process (i.e., 
typically on the night of observation or the following day). This is 
done using the GAMA-specific version of RUNZ, originally devel- 
oped by Will Sutherland for the 2dFGRS, and now maintained by 
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Figure 2. The flux correction vectors from each standard for one 2dF plate, 
shown here as an example. The grey lines show the individual flux correc- 
tion vectors for each standard star. The black line shows the B-spline fit to 
the mean vector. 



Scott Croom. This process is described in lDriver et al.l (|201 lh . and 
further details will be presented by Liske et al. (in prep), including 
the re-redshifting analysis, which quantifies the reliability of each 
measured redshift by having multiple team members re-measure 
redshifts for all low-quality flagged measurements, and for a subset 
of high-quality redshifts. 

The RUNZ code uses a cross-correlation approach to identify- 
ing an automated redshift, but allows the user to manually identify 
a redshift in the event of a poor result from the cross-correlation, 
before allocated a redshift quality flag (Q), from to 4, with 4 be- 
ing a certain redshift, 3 being probably correct, 2 indicating a pos- 
sible redshift needing independent confirmation, 1 indicating that 
no redshift could be identified, and meaning that the spectrum is 
somehow flawed and needs to be reobserved. Following this initial 
inspection, the process is repeated by multiple team members, in or- 
der to define a robust, probabilis tically defined "normalised" qual- 
ity scale, nQ jPriver et alj|201 1[). Detail s of th e error estimates on 
the redshifts are provided bv lDriver et al.ld201 lb . The re-redshifting 
process and derivation of the probabilities associated with the nQ 
quality will be presented in Liske et al, (in prep.). 

A small fraction (~ 3%) of the GAMA spectra are affected 
by fringing, as determined manually during the redshifting and re- 
redshifting processes. Of these ~ 50% still yield a good quality 
redshift, although other spectroscopic measurements such as emis- 
sion line properties (see §|5]below) are likely to be unreliable. 



4 SPECTROPHOTOMETRIC CALIBRATION AND 
QUALITY 

4.1 Flux calibration 

The main purpose of flux calibration is to first correct the 
wavelength-dependence of the system throughput (atmosphere, 
residual wavelength dependence of fibre entrance losses after the 
ADC, optics, and CCD quantum efficiency), and second to pro- 
vide an approximate absolute flux calibration. Obtaining accurate 
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Figure 3. (a) After dividing out the average flux correction vector (Fig. [2} 
to remove the low order shape, the low order residuals for each standard 
star are shown in grey. The residuals are then fit with a 4th order Legendre 
polynomial and the median of these polynomials is shown by the black line. 
The dotted line is the B-spline fit. (b) The high order terms are fit by first 
dividing out the low order (grey) and then taking the median of the result 
(black). 



spectrophotometry for the GAMA survey presents a challenge due 
to the 2" optical fibres used for spectroscopy, in addition to ob- 
serving in conditions that are not always photometric. Starting with 
the two-dimensional spectral output from 2DFDR, we spectropho- 
tometrically calibrate the data following the idlspec2d pipeline 
used for the SDSS DR6 jAdelman-McCarthv et al1l2008l) . We de- 
termine a curvature correction and relative flux calibration for each 
plate from the standard stars observed on each plate. The absolute 
spectrophotometric calibration is determined such that the flux of 
each object spectrum integrated over the SDSS filter curve matches 
the petrosian magnitude of the SDSS photometry for that object. 
We typically assign three of the fibres on each 2dF plate 
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Figure 4. The median ratio of common GAMA and SDSS spectra. This 
result is derived from 574 objects out of the 637 objects that were observed 
by both surveys. The 63 spectra excluded were due to particularly noisy 
GAMA spectra, or redshift mismatches between the surveys. The spectra 
are first normalised by the median flux value of the flux-calibrated spectrum, 
(since the absolute flux calibration is scaled differently between GAMA and 
SDSS). The spectra are then median filtered by 7 A and interpolated to the 
same wavelength scale before the ratio was taken. The dark line shows the 
median of the flux ratios, and the outer, grey lines show the 68th percentile 
range of the distribution of ratios for individual objects. 



for observing standard stars. The spectroscopic standards are typi- 
cally colour-selected to be F8 subdwarfs, similar in spectral type 
to the SDSS primary standard BD+17 4708. The spectrum of 
each standard star is spectrally typed by comparing it to a set 
of theoretical spectra generated from Kurucz model atmospheres 
(Kuru ca Il992l). using t he spectral synth esis code SPECTRUM 
I Gray & CorbaTivlll994l : iGrav et al.ll200lh . A llux correction vec- 
tor, a one-dimensional array of wavelength-dependent correction 
factors tied to the wavelength scale, is derived for each standard 
star by taking the ratio of its spectrum (in units of counts a nd af- 
ter correcting for Galactic reddening, ISchlegeletal1ll998h to its 
best-fit model (in units of ergss -1 cm -2 A -1 ). This is illustrated 
in Fig. |2] There are a small number of plates (2 out of a total of 
392 for v08 of the GAMA data) that included no standard star ob- 
servations. For these plates, we used the standard stars observed on 
the plate observed either just before or just after the plate lacking 
standard stars. 

We first calculate an average flux correction vector from the 
standard stars on each plate. We select high signal-to-noise regions 
of the standard star spectra, and divide out this average correction 
vector (Fig. [3^). By removing the overall average of the low-order 
shape in the standard star continuum in this way, we can fit the 
residuals to derive a plate-specific average curvature correction, to 
account for the declining CCD response at the extreme blue and 
red wavelengths. These residuals are fit with 4th order Legendre 
polynomials (Fig. [3^), and the average low order residual is found 
by taking the median of the Legendre coefficients. For spectral re- 
gions where the standard star spectra have S/N > 12, higher order 
fluctuations are also corrected. These higher order terms are found 
from the median Legendre coefficients after the low order terms are 
divided out. The flux calibration vector appropriate for each plate 
is constructed by multiplying the lower order residual fit and the 
high order residual fit. This vector is then fit with a B-spline, and 
the coefficients are used as the final flux calibration vector. Once 



the final flux calibration vector has been determined, it is applied to 
each individual spectrum on the plate, resulting in each spectrum 
being correctly flux calibrated, in a relative sense. 

The final step in obtaining an absolute flux calibration then 
involves tying the spectrophotometry directly to the r-band Pet- 
rosian magnitudes measured by the SDSS photometry. This is ac- 
complished by multiplying each individual spectrum by the SDSS 
r-band filter response. The SDSS magnitudes are based on pho- 
ton counting, so this calculation is done by integrating f v /v times 
the filter transmission function. We then determine the ratio of this 
flux with that corresponding to the SDSS Petrosian r-band magni- 
tude for each object. The GAMA spectra are then linearly scaled 
according to this value. It is straightforward to modify this scaling 
factor if alternative photometric references are preferred. In partic- 
ular, we are exploring the utility of Sersic r-band magnitudes as 
an alternative, although for the present the direct comparisons to 
SDSS measurements provide important consistency checks. 

To test how well the applied flux calibration method agrees 
with the methods applied to SDSS fibre spectroscopy, we look to 
objects observed by both surveys. There are 637 objects observed 
by both SDSS and GAMA, which are mostly galaxies, but also in- 
cluding 120 standard stars. After removing a small number of spec- 
tra to exclude particularly noisy spectra or those with mismatched 
redshifts, Fig.|4]shows the median and 68th percentiles of the ratio 
of the normalised GAMA and SDSS spectra as a function of wave- 
length, after the flux calibration is applied to the GAMA spectra. 
The spectra are normalised by the median flux of each spectrum 
before taking the ratio, due to the different absolute flux calibra- 
tions applied (to the fibre magnitude in SDSS, and to the Petrosian 
magnitude in GAMA). It is clear from the solid line that the flux 
calibration applied to the GAMA spectra results in a good agree- 
ment across the entire wavelength regime with high fidelity, and 
best in the mid-wavelength range of the spectra. The extreme blue 
and red ends of the spectrum are noisier, and worst in the blue, but 
the overall agreement resulting from the flux calibration applied to 
the AAOmega spectra from GAMA is still robust. Overall, we find 
an accuracy of around 10% in the flux calibrated spectra, ranging to 
somewhat worse than 20% at the extreme wavelength ends of the 
spectra. There remains an unresolved issue associated with the level 
of the response in the blue end of the spectra, which is the focus of 
ongoing work. For this reason, we advise caution when working 
with the bluest spectral diagnostics (such as the [Oil] emission line 
and the 4000 A break) in the current generation of measurements. 
Based on this analysis, we estimate that the flux calibrated GAMA 
spectra are typically accurate to between ~ 10 — 20%, although the 
small fraction of poorly spliced spectra and those affected by fring- 
ing are likely to be much worse. We note that this precision has 
been estimated using bright spectroscopic targets. For the fainter 
GAMA spectra it is likely that the spectrophotometric precision 
will not achieve this level, and we are also continuing to work on 
quantifying the dependence with magnitude. 

4.2 Spectral signal-to-noise 

Here we detail the continuum signal-to-noise (S/N) distribution 
for the GAMA spectra as a function of apparent magnitude (Fig. [5}. 
The S/N per pixel is quantified by measuring the median of the 
ratio of the observed flux to the noise in a 200 A window at the 
centre of the blue and red arms. The noise is determined from the 
variance arrays described in § 13.11 In the blue, the wavelength range 
4600 - 4800 A was selected, and 7200 - 7400 A in the red. It can 
be seen that the S/N properties vary as expected, with brighter tar- 
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Figure 5. The S/N in the blue (left) and red (right) arms of the GAMA spectra, as a function of g- or r-band Petrosian magnitude, respectively. The solid and 
dashed lines show the median and 68th percentiles. It is clear that brighter targets have higher S/N, and also that the red arm of the AAOmega spectrograph 
is relatively more sensitive for a given exposure time. The magnitude limit of r 17. 7 mag at the bright end corresponds to the SDSS main galaxy sample 
limit of r = 17.77 mag. The GAMA selection limit of r = 19.8 mag at the faint end is also apparent. The objects outside these limits i nclude bright target s 
not observed by SDSS due to fibre-collision limitations, and faint targets included through the additional K and 2-band selection limits (Baldry et al. 2010), 
together with specific additional filler targets (Ching et al. in prep.). 



Table 1. GAMA spectral emission lines measured by GANDALF. 



Emission line 


wavelength (A) 


Hell 


3203 


[NeV] 


3345 


[NeV] 


3425 


[Oil] 


3727 


[Nelll] 


3869 


H5 


3889 


[Nelll] 


3967 


He 


3970 


US 


4102 


H-y 


4340 


[oin] 


4363 


Hell 


4686 


[ArlV] 


4711 


[ArlV] 


4740 




4861 


[oin] 


4959 


[oin] 


5007 


[NI] 


5199 


Hel 


5876 


[OI] 


6300 


[OI] 


6364 


[Nil] 


6548 


Ha 


6563 


[Nil] 


6583 


[sn] 


6716 


im 


6731 



gets showing higher S/N. The red arm of the GAMA spectra dis- 
plays a typical S/N of a few at the faintest observed magnitudes, 
increasing to well over 10 at the brightest. In the blue arm the S/N 
is typically between 1 — 5 from the faintest to brightest sources. 
There are a very small fraction of spectra where the measured S/N 
in the blue is negative. This is a consequence of the continuum mea- 
surement being negative for these spectra, and is an artifact arising 
from poor flat-fielding or scattered light subtraction. There are also 



targets shown that extend fainter than our nominal survey selection 
limit of r = 19.8 mag. These enter the surve y as a result o f our sup- 
plementary K- and 2-band selection limits (Baldry et al. 2 010h . as 
well as the addition of a selection of "filler" targets. Due to galaxy 
clustering and the limitation of how closely fibres can be positioned 
with 2dF, there are frequently a number of fibres on any given plate 
that cannot be assigned to a primary survey target. In this case, and 
in order to maximise the scientific return from the survey, we allow 
fibres to be assigned to supp lementary targets selected to support 
specific scientific programs dDriver et al .11201 ll) . In particular, we 
allocate such filler fibres to target systems identified as interesting 
objects in the far-infrared by the Herschel-KILAS survey, as well 
as to radio sources with carefully identified optical counterparts, 
that otherwise lack redshift or spectroscopic measurements (Ching 
et al. in prep.). 



5 MEASUREMENT OF SPECTRAL FEATURES 

The flux-calibrated ID spectra are those used for all subsequent 
measured spectroscopic properties. Both emission lines and ab- 
sorption features are measured. A variety of standard galaxy prop- 
erties are catalogued based on these measurements, including 
B aimer decrements and star formation rates from the Ha and 
H/3 emission lines i Gunawardhana et al]|2013l) . spectral diagnos- 
tics jBaldwin et alj|l98ll) to discriminate AGN from star forming 
objects dGunawardhana et al ] |2013l) . metallici y measurements for 
the nebular gas from the emission line features dFoster et al J|20 1 2t 
Lara-Lopez et al., in prep.), stellar velocity dispersions from ab- 
sorption features and the D4000 age-estimate parameter from the 
4000 A break, and more. This section details the measurement pro- 
cesses for quantifying the spectral features used in deriving these 
now standard properties for GAMA galaxies. The measurement of 
the specific properties themselves are detailed in the various papers 
presenting the analysis of those parameters, referenced above. 

Emission lines are measured in two ways. We first fit Gaus- 
sians to a selection of common emission lines at appropriate ob- 
served wavelengths, given the measured redshift of each object. 
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Figure 6. Upper left: The distribution of the ratio of [OIII]A4956/[OIII]A5007, showing consistency with the value fixed, for a giv en density and temperatu re, 
by quantum mechanics. The solid line is the median, and the dashed lines are the 68th percentiles. The expected ratio of 1/2.98 jStorev & Zeippenll2000l) is 
shown as the dot-dashed line. Upper right: The flux ratio shown now as a function of the flux in the brighter line. The dashed line here is the expected value. 
The bottom panels reproduce this analysis for [SII] A6716/[SII]A6731, where the expected ratio is 1/1.4 (bsterbrock 1989b . 



This is performed for a subset of common emission lines via a si- 
multaneous iterative \ 2 fitting of positive emission peaks to these 
emission lines in three independent line groups, around the wave- 
lengths of the [Oil], H/3, and Ha lines. The local continuum span- 
ning each fitting region is approximated with a linear fit. All lines 
in each group are fitted simultaneously. As well as the continuum 
coefficients and peak intensities, a small velocity offset (from the 
underlying GAMA redshift) is allowed, and a line width common 
to all lines in the group is also fit. Line width is constrained to 
lie between 0.5-5 times the instrumental PSF (3.4 CCD pixels, 
3-5A). Each line group is fitted with independent values for the 
free parameters to accommodate small local variations in the spec- 
tral sampling. Flux values for individual lines are rejected if in- 
clusion of a line in the global fit fails to improve the reduced \ 2 
value significantly (by a factor of 3). Typical GAMA sources ex- 
hibit only marginally resolved emission lines at the resolution of 
the AAOmega spectra. No attempt has been made to fit multiple 
emission components to composite emission line structures, such 
as the [Oil] doublet, at this time. Integrated line flux, EW and the 
associated statistical errors are estimated from the fitting process. 
The errors are those associated with the formal Gaussian fitting pro- 
cess, and as shown in § 16. 21 below, the error estimates are robust 
although likely to be somewhat underestimated for fainter galax- 
ies. In addition, flags are provided to identify a variety of problem- 



atic cases, including failed measurements or lines falling in masked 
spectral regions, lines that are apparently too narrow to be fit (usu- 
ally caused by bad pixels), and lines that are fit by the maximum al- 
lowed width of the fitting routine (a Gaussian a — 10 A), typically 
due to the presence of an intrinsically broad line. This provides a 
baseline set of emission line fluxes, equivalent widths, and S/N 
estimates. For those (small fraction of) objects such as broad-line 
active galactic nuclei (AGN) and sources that exhibit line splitting 
(indicative of either merging systems or starburst winds), the sim- 
ple single Gaussian emission approximation will not provide useful 
estimates. 

Independently, we also pursue a more sophisticated spectral 
measurement approach, in order to self-consistently derive both 
stellar kinematics and emission line properties. To do this, we use 
the publicly avail able codes pPXF dCappellari & Emsellemll2004h 
and GANDALF JSarzi et al.l2006l) . We first extract the stellar kine- 
matics using pPXF, matchi ng the observed spectra t o a se t of stellar 
population templates from Maraston & Stromback (201 1) based on 
the MILES stellar library dSanchez-Biazquez et al.l2006h . masking 
the regions which could potentially be affected by nebular emis- 
sion, given the observed redshift. For this step, we dow ngrade the 
spectral resolution of the lMaraston & S tromback ( 1201 lh models to 
match the GAMA spectral resolutio n, adopting a value of 2.54 A 
(FWHM) for the MILES resolution dBeifiori et alj201lh . Next, we 
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use GANDALF vl.5 to simultaneously fit both Gaussian emis- 
sion line templates and the stellar population templates, which are 
broadened to account for the derived stellar kinematics, to the data, 
while also correcting for diffuse (stellar continuum) obscuration. It 
is important to note that diffuse obscuration, i.e. that caused by dif- 
fuse dust in the galaxy, affects the entire spectrum, and no obscura- 
tion correction is applied purely to emission lines during the fitting 
process (even though this option is available). We made this de- 
cision in order to minimise the number of poor spectral fits caused 
by low spectral S/N that otherwise led to ambiguous emission-line 
fluxes, when n oise was a more significant factor than reddening. A 
ICalzettild200lh obscuration curve is used in estimating the obscura- 
tion corrections to the stellar continuum, and is applied throughout 
the GANDALF measurement process. The nebular emission lines 
from GANDALF are subsequentl y corrected for obscu ration using 
a Milky W ay obscuration curve dCardelli et al.l ll989). as recom- 
mended bv lCalzettilj200lh . 

GANDALF's simultaneous fitting mechanism allows us to ac- 
curately extract ionised gas emission from the stellar continuum, 
minimising contamination from stellar absorption, in order to cal- 
culate emission line fluxes and gas kinematics (velocities and ve- 
locity dispersions) from the Gaussian emission line templates. The 
outputs from this analysis are line fluxes and equivalent widths, 
velocity dispersions (from the stellar absorption lines in the best- 
fitting SEDs), and associated derived products such as Balmer 
decrement. A list of emission lines for which flux and equivalent 
width measurements are extracted is given in Table [T] The wave- 
lengths given for these lines are the wavelengths as meas ured in air, 
as op posed to the vaccuum wavelengths used by SDSS dYork et all 
l200d) . 

Our final set of measurements includes the flux, equivalent 
width and signal-to-noise ratio for each emission line, as well as 
a velocity dispersion inferred from the line w idth. We calculate an 
emission line ratio diagnostic classification jBaldwin etai]|l98lh 
for each galaxy based on these measurements. The stellar veloc- 
ity dispersion and E(B — V) values from the SED fits, for both 
diffuse and nebular obscuration (if applicable), are also recorded. 
While stellar velocity dispersions are measured for all spectra, 
these are typically only robust for spectra having high S/N (e.g., 
IProctor et al.ll2008l : IShu et al.ll201 ll Thomas et al., in prep.). In ad- 
dition, the best-fit SED template and associated y/ 2 > along with a 
clean emission-line free absorption spectrum are available. 



6 QUANTIFYING MEASUREMENT RELIABILITY 

We explore the quality of our spectroscopic measurements thor- 
oughly. Here we detail the internal reliability for measurements 
within each spectrum (§ 16. lb . the repeatability of our measurements 
using duplicate observations (§ 16, 2b . and the self-consistency of our 
measurements between the two independent approaches we use 
(§ 16. 3b . We go on to outline the impact of stellar absorption on the 
Gaussian fit Balmer line measurements (§ !6.4b , the reliability of ve- 
locity dispersion estimates (§ !6.5b . and the extent of aperture effects 
(MM- 



6.1 Internal consistency 

A simple test of the internal consistency of our line measurements 
is to measure the ratios of emission lines from ionised species 
that should be fixed by quantum mechanics for fixed density and 



temperature, with the [OIII], [Nil] and [SII] line pairs being ob- 
vious choices. Due to the stellar absorption of Ha in the vicin- 
ity of the [NII]A6548 line, and (to a lesser degree) its proxim- 
ity to the Ha emission line itself, this particular ratio is less ro- 
bust. Using the measurements from the Gaussian fits, in Fig. [6] 
we show the distribution of [O III] A4956/[OIII] A5007 c ompared to 
the expected ratio of 1/2.98 dStorev & Zeip pen 2000). The same 
is sho wn for [SH]A67 16/[SII]A6731, where the expected ratio is 
1/1.4 JOsterbrockll 1989b . Line ratios are only included in this anal- 
ysis for pairs where the brighter of the two lines has a flux mea- 
surement above 3cr. In the right panels, the ratio is shown as a 
function of the flux of the brighter line. At lower flux (more specif- 
ically lower S/N), the fainter line is less easily detected, and a bias 
toward higher values of the ratio can be seen. This result demon- 
strates the robustness of the line fitting and measurement within 
each spectrum. 

6.2 Duplicate measurements 

A small number of GAMA spectra duplicate observations of partic- 
ular targets, often due to the re-observation of objects where a low 
redshift quality was initially obtained. This is typically a conse- 
quence of the objects being at the fainter end of the GAMA target 
selection, and the results here are consequently illustrative of the 
robustness of the measurements for the fainter population. There 
is no attempt made to combine these duplicate GAMA spectra, or 
the measurements from them, due to the complex systematics in- 
volved in the flux calibration steps detailed above. We can, how- 
ever, take advantage of these duplicate observations to understand 
the precision to which we can measure our emission line equiva- 
lent widths and fluxes. To do this, we select those duplicate mea- 
surements wh ere both spectra hav e been allocated a redshift quality 
3 ^ nQ ^ 4 toriver et al.ll201 ll) . We illustrate the differences be- 
tween these duplicate measurements, using both the GANDALF 
and the Gaussian fits, for line fluxes (Fig. [7} and equivalent widths 
(Fig.(8]( of the Ha and [Oil] emission lines. These two line species 
were chosen for the following reasons. First, they represent the 
best- and worst-case scenarios with Ha being (typically) a strong 
line in a high S/N region of the spectrum, and [Oil] being both 
weak (in many cases) and situated at the low S/N end of the spec- 
trum. Thus, these two lines give an idea of the spread expected 
in the precision of the measurements across the full wavelength 
and S/N range. Second, the Ha line is affected by underlying ab- 
sorption due to the stellar continuum which is corrected for dur- 
ing the GANDALF fitting process. However, this correction relies 
critically on a robust model fit to the underlying stellar continuum 
and errors in this fit may add systematic uncertainties to the flux 
and equivalent width measurements. Such systematics are quanti- 
fied by comparing the differences in the duplicate measurements 
to the quadrature sum of the uncertainties on the measurements 
(Fig. [9] and Fig. 1 10b. The 68th percentiles are typically 1.5 to 2 
units from the median, suggesting that the formal statistical errors 
on the line measurements are somewhat underestimated from both 
the GANDALF and the Gaussian fits. This is likely related to the 
fact that the duplicate spectra are dominated by fainter galaxy tar- 
gets. For the fainter or lower S/N targets, the continuum level in 
the spectra tends to be noisier, and more affected by systematics 
such as poor scattered light subtraction or sky subtraction. This 
leads in turn to greater variation between the repeat measurements, 
as the continuum level estimated for either the Gaussian fitting or 
by GANDALF can be more easily over- or under-estimated. The er- 
rors on the line fitting may not accurately reflect the uncertainty in 
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Figure 7. Top row: The distribution of the ratio of the differences in the GANDALF emission line fluxes to the mean of the two flux measures for [Oil] (left) 
and Ho (right), for repeat measurements of the same galaxy. The difference A/ = /i — /a is taken in the sense that fi comes from the spectrum with the 
higher S/N of the pair of flux measurements. Bottom row: Distribution of the relative differences for Gaussian fit measurements. In each panel the solid line 
shows the median of the distribution, and the dashed lines show the 68th percentiles. 
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Figure 8. Same as Fig.|7]but for [Oil] (left) and Ha (right) equivalent widths rather than fluxes. The top row shows the GANDALF results, while the bottom 
row shows the Gaussian Fit results. 
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Figure 9. Same as Fig.[7]but the difference is normalised by ei,2 = (<rj + If no biases in the measurements exist due to differences in the spectrum 
S/N ratio, then the distributions should be centred at zero and be symmetric. If the uncertainty estimates are accurate measures of the true uncertainties, then 
the 68th percentile values should be of order unity. 
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Figure 10. Same as Fig.|9]but for the line equivalent width measurements, rather than the fluxes. 
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Figure 11. Comparison of line flux measurements using GANDALF and Gaussian-fit methods. The top row shows the direct comparison of the two methods, 
with the deviation from the one-to-one relation for the Balmer lines being due to stellar absorption. The bottom row shows a histogram of the differences 
in the measurements normalised by the GANDALF measurement. The tails illustrate cases where the two values are quite different, with the spike visible 
at unity corresponding to cases where the Gaussian-fit flux is significantly smaller than the GANDALF value. The tail to negative values corresponds to the 
GANDALF measurement being significantly smaller than the Gaussian-fit value. 
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Figure 12. Same as Fig. 1121 but for line equivalent widths rather than fluxes. 
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the continuum estimation. For brighter (higher S/N) targets, this is 
less likely to be a limitation for the Gaussian fits, although system- 
atic underestimates in the uncertainties may still be possible in the 
GANDALF fitting if the continuum is not well-described by the 
underlying SEDs. Overall, the repeatability of the duplicate mea- 
surements is very high, with 1 a differences of less than 0.5 A for 
the equivalent width of Ha and less than 1 A for [Oil]. The re- 
peatability of the flux measurements is similarly reliable, with 1 a 
variations of typically less than 0.5 x 10 -17 ergs -1 cm -2 . 



6.3 Self-consistency 

Having established the precision of the emission line measurements 
in the GANDALF and Gaussian-fit catalogues individually, us- 
ing the duplicate measurements, we now test for systematic biases 
which may be inherent in the different techniques used to determine 
the line equivalent widths and fluxes. To identify such potential bi- 
ases, in Fig. QT] and Fig. [T2] we investigate comparisons between 
the two measurement methods for the equivalent widths and fluxes, 
respectively, for the [Oil], H/3, [OIII] A5007, Ha and [Nil] A6584 
line species. In these Figures, we only include comparisons where 
the emission line of interest is detected with an amplitude-to-noise 
ratio A/rN > 3 as determined for the GANDALF fits. Here 
A/rN is the ratio of the line amplitude (from t he Gaussian fit) t o 
the standard deviation of the residual spectrum (ISarzi e t al. 2006). 
We also limit ourselves to cases where both line measurements have 
an equivalent width greater than zero (noting that we take the con- 
vention that emission lines have positive equivalent widths). The 
forbidden lines, [Oil], [OIII] A5007 and [Nil] A6584 are chosen 
because they are not significantly affected by stellar absorption and 
probe a significant portion of the wavelength range covered by the 
AAOmega spectra. Since the GANDALF measurements correct for 
the effects of stellar absorption on the emission line measurements, 
the comparison between the Balmer line species, H/3 and Ha, al- 
lows us to quantify the systematic effects of the underlying stellar 
absorption on the Gaussian-fit measurements. 

Considering the forbidden lines only, the distributions are con- 
sistent with the one-to-one relation. The median and modes of the 
relative difference in the fluxes for [Oil] and [OIII] A5007 (lower 
panels, Fig. II It indicate offsets less than 5% (with a slight system- 
atic toward the GANDALF measurements being larger). For the 
equivalent widths of [Oil] and [OIII] A5007, however, the median 
and modes of the relative differences (lower panels, Fig. I12t indi- 
cate that the GANDALF values are systematically higher by around 
10 — 15%. The likely cause of this small offset in the equivalent 
width measurements is in the different definitions used for the con- 
tinuum flux estimates. The [Nil] line, though, does not appear to 
be affected by an offset of the same magnitude (equivalent width 
differences of at most 5%), which further suggests that the bulk of 
the difference arises in the continuum estimate for the noisier blue 
arm of the spectra ([Oil] is always in the blue for GAMA spectra, 
while [OIII] is in the blue for spectra with z < 0.13). The con- 
clusion here is that the independent measurement of emission line 
fluxes are consistent in the median to better than 5%, with a disper- 
sion consistent with the error measurements on the lines. The line 
EW estimates are also consistent to better than 5% in the red, and 
to 10 — 15% in the blue, again with dispersions consistent with the 
measured errors. 
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Figure 13. Balmer decrement (Fua/FHp) as measured by GANDALF, 
compared against that from the Gaussian line fits after correcting for a con- 
stant stellar absorption equivalent width of 2.5 A. The solid line indicates 
equality, and the dashed lines correspond to the value of 2.86 expected for 
Case B recombination. 



6.4 Stellar absorption 

As expected for the Balmer lines, the GANDALF flux and equiva- 
lent width measurements are systematically higher than the Gaus- 
sian fit measurements due to stellar absorption, which is accounted 
for by GANDALF but not in the Gaussian fits. The difference is 
particularly conspicuous for the weaker H/3 emission line where 
the stellar absorption can be large compared to the typical H/3 
emission strength. The relative difference in the equivalent width 
measurements for H/3 is offset by ~ 50%, while the offset is 
~ 15% for Ha. The offset is dominated by the correction to the 
underlying stellar absorption in the GANDALF measurements, al- 
though the systematic offsets affecting the equivalent width mea- 
surements for [Oil] and [OIII] may also be present at some level. 
We find that applying an average stellar absorption correction to 
the equivalent widths of 2.5 A is appropriate in order to make 
the Gaussian line fits consistent with those from the GANDALF 
measurements. This can be seen explicitly in a comparison of the 
Balmer decrements from the GANDALF measurements compared 
against the stellar-absorption-corrected Gaussian fits in Fig. [13] 
This figure also indicates the typical Case B recombin ation value 
of 2.86 for the Balmer decrement dOsterbrockl Il989l) . assuming 
T e = 10000 K and n e = 100 cm" 2 . Note, though, that the in- 
trinsic Balmer decrement c an be as high as ~ 3 for temper atures 
T e 5000 - 6000 K (e.g jLopez-Sanchez & Estebanll2009l) . The 
inferred stellar absorption correction of 2.5 A is somewhat smaller 
(by ~ 1 A) than the typical stell ar absorption equivalent widths 
for disk galaxies (Kenn icutill 19831) . but consistent with the result 
found by Hop kins et alj J2003h 7 This is also consistent with stellar 
absorption equivalent widt hs, from 1 — 2.5 A, found in a s ample 
of star forming galaxies by Lopez-Sanchez & Esteban (20 10[), who 
show that the stellar absorption tends to increase with increasing 
stellar mass, or metallicity, of a galaxy. When the stellar absorption 
is marginally resolve d, as it was for the SDSS spectra analysed by 
Hop kins etall feoOS). and as it is for the current AAOmega spectra, 
a Gaussian fit to the line flux requires a smaller correction. 
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Figure 14. Differences between repeat measurements of the velocity dis- 
persions, normalised by the quadrature sum of the errors. This analysis is 
restricted to duplicate spectra where both of the pair have S/N > 3 in the 
blue. The 68th percentile of this distribution is about ±2, suggesting that 
the errors in the velocity dispersions in these relatively low S/N spectra 
are underestimated by a factor of two. 




Figure 15. Effective aperture correction as a function of redshift. The cor- 
rection is a multiplicative flux scaling. The colour coding indicates the stel- 
lar masses of the galaxies. The majority of GAMA systems have aperture 
corrections of 2 — 4, with a small number of galaxies at the lowest redshifts 
having larger corrections. 



6.5 Velocity dispersions 

Velocity dispersions, as measured by pPXF, are reliable for rel- 
atively high signal-to-noise spectra. Earlie r work suggests a con- 
servatively high threshold of S/N > 12 jProctor et al1l2008l) . for 
extracting reliable velocity dispe rsions, but more recent work sug- 
gests that spectra with S/N > 5 dShu et alj201 ll Thomas et al., in 
prep.) may still yield reliable measurements. Fig. [l4]shows the re- 
peatability of the velocity dispersion measurements from duplicate 
observations. The 68th percentile of this distribution is around two, 
implying that the measured velocity dispersion errors are underes- 
timated by about a factor of two. This is likely a consequence of 
the relatively low S/N of the duplicate spectra available to make 
this measurement, though, especially in the blue (Fig. [5}, since the 
absorption features used in constraining the velocity dispersion lie 
primarily in the blue half of the spectrum. In particular, errors in 
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Figure 16. Balmer decrement as a function of aperture corrected (but not 
obscuration corrected) Ha luminosity. Systems classified as AGN are not 
shown (see § |7.2| below). The dashed line shows the Case B recombination 
value of 2.86, and the solid line shows a fit to the observations with Balmer 
decrements higher than this value. 



the estimate of the continuum level are again likely to be contribut- 
ing substantially to the underestimate in the errors on the velocity 
dispersions, for the low S/N spectra. We have used duplicate spec- 
tra with S/N > 3 in order to sample enough duplicates for this 
analysis, since limiting the analysis only to S/N > 12 rejects the 
majority of the duplicate spectra. Consequently, since the veloc- 
ity dispersion measurements rely largely on absorption lines at the 
blue end of the spectra, these low S/N duplicate spectra are likely 
to be less well characterised than those of higher S/N. 



6.6 Aperture effects 

Due to the 2" fibre diameter, emission from a galaxy that is larger 
than this on t he sky will not be measured within the fibre (e.g., 
Hopkins et al. 2003; Brinch mann et afl l2004l) . A detailed analysis 
of the systematic effects in fibre-spectrograph measurements asso- 
ciated with fibre-po sitioning e rrors, efficiencies and aperture cor- 
rections is given by|N ewm an] J2002I) . Below we outline the extent 
of the aperture effects in the GAMA spectra. 

Corrections to account for these aperture losses are incor- 
porated into the flux calibration process above (§0, and are ef- 
fectively multiplicative flux scalings. These can only ever be ap- 
proximations, of course, making the assumption that the spectro- 
scopic properties (such as line fluxes, or derived properties such 
as star formation rate) that lie outside the fibre aperture can be 
scaled by the broadb and light profile available from the photometry. 
iKewlev etalJ (2005) highlight, for example, how aperture effects 
can bias the estimate of star formation r ate, n ebular metallicity and 
obscuration in galaxies. iGerssen et alj (120121) presents an analysis 
highlighting the potential biases involved, for a sample of galax- 
ies at z < 0.1, in particular comparing estimates of star formation 
rate based on integral field data to those from aperture-corrected 
single-fibre measurements. They emphasise that a simple aperture 
correction of the kind applied here can underestimate the true ef- 
fect by factors of ~ 2.5 on average, although with a large scatter. 
Circumventing the limitations of such a fixed fibre covering frac- 
tion can only truly be achieved with multiobject integral field spec- 
troscopy, pro mised by the next g eneration of instrumentation such 
as the SAMI dCroom e t al. 2012) instrument on the AAT. There is 
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Figure 17. Balmer decrement as a function of aperture corrected (but not obscuration corrected) Ha luminosity. The top row shows the diagram separated into 
bins of redshift, the bottom row into bins of galaxy stellar mass, both increasing left to right. The dashed line is again the Case B value of 2.86. 



no simple systematic solution, because of the underlying variety of 
differently distributed star formation locations within galaxies. In 
the absence of such detailed measurements, a simple aperture cor- 
rection still remains the best proxy for a total line luminosity or star 
formation rate for a galaxy, and we briefly discuss here the scale of 
these corrections for the GAMA spectra. 

To illustrate the extent of the aperture corrections, Fig. [T3] 
shows the aperture correction as a function of redshift. This can 
be compared against a perture corrections u sing a similar approach, 
seen in Fig. 25(a) of Hop kins et alj d2003l) . which shows the same 
quantity (given logarithmically) as a function of redshift for SDSS 
spectra. For the SDSS galaxies, the aperture corrections vary from 
factors of ~ 2 to 10 or more at the lowest redshifts. Comparing 
the SDSS galaxies with the GAMA targets, a larger fraction of the 
SDSS systems below z ~ 0.1 have aperture corrections larger than 
a factor of 5, due to the larger size on the sky of the low-redshift 
galaxies, despite the larger SDSS fibre aperture. 

The GAMA spectra have been taken using 2" diameter fi- 
bres, compared to the 3" diameter fibres used by SDSS. This is 
not as problematic as might be initially assumed, however, as the 
GAMA targets (primarily 17.77 < r < 19.8) are fainter than those 
from SDSS (r < 17.77). Consequently, at a fixed redshift the typ- 
ical GAMA target is significantly smaller (in kpc) than the typical 
SDSS galaxy. Similarly, at fixed mass, the typical GAMA target 
is at higher redshift, and so smaller (in arcsec). It turns out that the 
distribution of aperture sizes in units of the effective radius are very 
similar between GAMA and SDSS. Thus, as a consequence of the 
fainter magnitude limit, and greater redshift depth, despite the fact 
that 2dF/AAOmega has smaller fibres than SDSS, we are no more 
susceptible to aperture effects. 



7 EMISSION LINE RATIOS 

Ratios of bright emission line measurements for galaxies are com- 
monly used to constrain the obscuration (through the Balmer decre- 
ment, Hq/H/3), or to discriminate between a supermassive black 
hole (active galactic nucleus, AGN) or star formation (SF) as the 
photoionisation source. A typical spectroscopic diagnostic use for 
the latter discrimination compa res the [OlIl]/H/3 to [Nii]/Hq ratios, 
following lBaldwin et alj 1 1198 II) . In this section we demonstrate the 
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Figure 18. Spectral diagnostic diagram demonstrating the range of diagnos- 
tic m easurements in the G AMA sample. Th e discrimination lin es shown are 
from lKewlev et alj (200 II) (solid) and Kauffmann et al. 1 2003b) (dashed). 



range of such properties present in the GAMA sample. Here we 
show results using the Gaussian fits to the emission lines, although 
we see identical trends if we use the GANDALF measurements. 



7.1 Balmer decrement 

The Balmer decrement is calculated as the ratio of stellar- 
absorption corrected line flux, Ha/H/3. The distribution of Balmer 
decrement, as a function of stellar absorption corrected and aper- 
ture corrected Ha lumin osity, is shown in Fig. [16] (see also 
iGunawardhana et all 1201 3l) . Galaxies are only shown in this fig- 
ure if they have Sua > 25 x 10 -17 erg s" 1 cm~ 2 , and local 
volume flow corrected redshifts in the range 0.001 < z < 
0.35. The flux limit corresponds to the limit at which the emis- 
sion lines are robustly measured jGunawardhana et al.ll2013l) . The 
lower redshift limit excludes galaxies with e rroneously low red - 
shift measurements or stellar contaminants fealdrv et alj l2012h . 
The upper redshift limit is where Ha falls outside the observ- 
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able spectral range. The trend demonstrates that high line lumi- 
nosity systems show a much broader distribution of obscuration 
properties than lower luminosity ga l axies, consistent with ear- 
lier results (H opkins e t al.l200lLl2003l:|Perez-Gonzalez et alj2003l; 
Afonso et alj 120031; iLopez-Sanchezj l20ld: iGunawardhana et al.l 



201 lh . 



With the numbers of galaxies available to us in the GAMA 
sample, we can explore the distribution of obscuration properties 
in more depth by looking at trends with both mass and redshift. We 
find results that are consistent with, and c omplementary to the early 
SDSS work bv lKauffmann et alj <2003al l as well as many other re- 
cent works. Fig.QTJdemonstrates that the low stellar mass galaxies 
tend typically to have lower Ha luminosities, and low overall lev- 
els of obscuration. Progressively higher mass systems, which can 
sustain higher levels of star formation, tend to display both higher 
levels of Ha luminosity and more extreme Balmer decrements. 
Note that Fig. [T7] shows the Ha luminosity before correcting for 
obscuration, so the difference in intrinsic luminosities between low 
and high mass systems will be enhanced. Interestingly, at the high- 
est masses, there is a detectable division into high Ha luminosity, 
high obscuration systems, and lower Ha luminosity, low obscura- 
tion systems. These high-mass, lower luminosity systems have low 
specific star formation rates, and given their low obscurations are 
likely to correspond to objects underg oing the transitio n from the 
blue cloud to the red sequence (e.g. jBaldrv et al.ll2004l) . This pop- 
ulation is explored in more detail in Taylor et al. (in prep). 

When these trends are explored as a function of redshift the 
same broad picture emerges, a consequence of the apparent magni- 
tude limit of GAMA resulting in high-mass systems being prefer- 
entially identified at higher redshift. There is a detectable difference 
at higher redshift, with the absence of high-mass systems display- 
ing low luminosities and obscurations. This is likely to be a con- 
sequence of Malmquist bias, with such low-luminosity systems not 
being able to be detected at higher redshift. 

The interesting point to highlight here is that there is a bivari- 
ate selection effect at work when exploring the properties of emis- 
sion lines in a spectroscopic survey of a broad-band magnitude- 
limited sample. The broad-band magnitude limit, to first order, cor- 
responds to a (redshift-dependent) stellar-mass limit. This is subse- 
quently subjected to an emission line flux limit through the spec- 
troscopic observations, which, again to first order, corresponds to a 
(redshift-dependent) luminosity limit. There will always be galax- 
ies that are bright enough in the continuum to be targeted spectro- 
scopically, but whose emission properties are too faint to detect. 
These can be accounted for with appropriate completeness correc- 
tions. There will also exist galaxies that are brighter than the spec- 
troscopic sensitivity limits, but which never enter the spectroscopic 
sample as their hosts are fainter than the broad-band magnitude lim- 
its used to select the targets. This component cannot be accounted 
for with completeness corrections since their host population is not 
well-defined. The consequences of these bivariate selection effe cts 
in GAMA are discussed in detail in lGunaw ardhana e t al.l ( l2013l) . 



7.2 Spectral diagnostics 

The ratios of forbidden emission lines to Balmer lines have 
been used for many decades now as discriminators be- 
tween different sources of photoionisation ( Baldwin et al.l [l98ll : 



10~ 17 ergs~ 1 cm" 2 , and local volume flow corrected redshifts in 
the range 0.001 < z < 0.35. This Figure discriminates star form- 
ing galaxies (below the dashed line) from those where the ioni- 
sation arises from an AGN (above the solid line), with galaxies 
between these discriminators commonly treated as composite sys- 
tems. Of the emission line systems in GAMA, the majority are 
classified as star forming in this fashion. Only about 12% of the 
galaxies with measured [Nil], [OIII], Ha and H/3 fluxes, and qual- 
ity nQ ^ 3, are classified as AGN. 

Again capitalising on the sample size available with GAMA, 
we explore the mass and redshift dependencies of the star forma- 
tion and AGN distributions, shown in Fig. [19] Perhaps not too sur- 
prisingly, the trends visible here highlight that AGN systems are 
more prevalent in more massive galaxies, and given the magnitude- 
limited nature of the sample, these are more visible at higher red- 
shift. At the lowest masses, the galaxy population is almost com- 
pletely dominated by star forming systems. As galaxy stellar mass 
increases, there is a progressive increase in the proportion of AGN 
systems. T hese results are consistent with those dem onstrated from 
the SDSS jKauffmann et al]|2003b! : lHao et alj2005l) . 

It is also worth noting the morphology of the region populated 
by the star forming galaxies in Fig. [19] The star forming population 
moves from a region of low [NII]/Ha and high [OIIIj/H/? for low 
mass systems, corresponding to low metallicities, progressively to 
having high values of [NII]/Ha and low values of [OiIIj/H/3 for 
galaxies of high stellar mass, corresponding to high metallicities. 
This transition reflects th e well-established mass-metallicity rela- 
tionship for galaxies (e.g-lTremonti et al.ll2004l ; lKewlev & Ellison! 
I2008L lLara-Lopez et al.l l201dh . Details of the metalli city proper- 
ties of gala xies in the GAMA sam ple are presented bv lFoster et al.l 
j2012l) and lLara-Lopez et al.l J2013I) . 

The redshift dependencies of this spectral diagnostic are also 
illuminating. As with the Balmer decrements, in a broad sense the 
redshift trends reflect the mass dependencies due to the magnitude- 
limited nature of the survey, which leads to high-mass systems 
preferentially being found at higher redshift. Interestingly, though, 
at the lowest redshift there is a populati on of LINER-like AGN , 
with high [NII]/Ha with low [OIII]/H/3 dSchawinski e"tafll2007l) . 
which are relatively high-mass systems. This ionisation signa- 
ture may also be more characteristic of shock-excitation than the 
mo re prevalent active nucleus driven ionisation in massive galax- 
ies JSharp & Bland-Hawthornll2010l : iFarage et al.ll2010b iRich et al.l 
|2010|.|2011|) .~ 



_ _ l_pl 

IVeilleux & Osterbrockl fl987t iKewlev et al]|200lh . In Fig. [T8l we 
show the most commonly used diagnostic diagram, [OIII]/H/3 
as a function of [NII]/Ha, for the full GAMA sample. Sys- 
tems shown here are again limited to those with Sua > 25 x 



8 CONCLUSION 

We have detailed above the processes involved in compiling, pro- 
cessing, calibrating and measuring the AAOmega spectra that un- 
derpin the GAMA multiwavelength survey. Details of the spectro- 
scopic flux calibration provided in §[4]show that we achieve a preci- 
sion of about 10—20%. The continuum S/N in the spectra is higher 
in the red than in the blue, being typically ~ 10 in the red and ~ 5 
in the blue for the brightest targets, and decreasing as expected for 
fainter targets. The spectroscopic measurement reliability has been 
quantified in terms of internal consistency, repeatability and self- 
consistency between independent approaches to the emission line 
measurements. These analyses demonstrate that the various mea- 
surements give consistent results with robustly estimated uncertain- 
ties. It is important to note that we provide both the relatively di- 
rect Gaussian fit measurements, as well as those from the more so- 
phisticated GANDALF fitting, as GAMA data products. We do this 
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recognising that there will be some spectra for which GANDALF 
is not able to make a reliable measurement (such as a badly spliced 
spectrum where the SED fitting has failed), but for which reliable 
Gaussian fits to the emission features can still be made. There are 
also likely to be a variety of science cases where having a larger 
number of simple measurements is more valuable than having a 
smaller number of more refined measurements, and vice-versa. To 
facilitate both aspects, both sets of measurements are provided. 

The GAMA survey has already produced a broad cross- 
section of insights into the properties of galaxy evolution, and as the 
survey progresses it will continue to provide a unique and valuable 
spectroscopic and multiwavelength resource for studies of galaxy 
formation and evolution for many years to come. All the GAMA re- 
sults that have been published to date are based only on data from 
GAMA I, i.e., the 144 deg 2 contained within the three Equatorial 
fields, G09, G12, G15, taken during observing campaigns spanning 
2008-2010. The GAMA survey has continued through 2010-2012 
with additional observations of the Equatorial fields, to expand the 
area and achieve a uniform survey depth of r pe t < 19.8 mag, as 
well as opening up two new Southern fields, G02 and G23. The 
goal is to survey a total area of 280 deg 2 to a uniform depth of 
?"pct < 19.8 mag, resulting in ~ 300 000 galaxy spectra. To date, 
over 220 000 spectra have been measured. 

We illustrate the data quality and utility with a simple ex- 
ploration of how obscuration in galaxies varies with galaxy mass 
and redshift, using the Balmer decrement. We find, consistent with 
earlier work, a luminosity-dependence in galaxy obscuration. This 
effect is seen both as a function of mass and redshift, largely as 
a consequence of the magnitude-limited nature of the survey. We 
do identify, though, a population of high mass, low Ha luminos- 
ity systems, with relatively low obscuration, that are likely to be 
systems transitioning from the blue cloud to the red sequence. We 
also explore the mass and redshift dependence of the spectral di- 
agnostic diagram, finding that AGN systems are more prevalent 
in higher mass galaxies, which are more numerous at higher red- 
shifts in magnitude limited samples like GAMA. Higher mass sys- 
tems are also less likely to have lower [NII]/Hq ratios, consistent 
with having higher nebular metallicities, and reflecting the well- 
established mass-metallicity relationship for galaxies. Both of these 



results and more are being explored in more detail in a variety of 
works in progress. 

The raw and processed GAMA spectra, and the de- 
rived data products, are being progressively released to 
the public through staged Data Releases. The data and 
data products will be available from the GAMA web site 
|http : //www . gama-survey . org/| The GAMA team 
welcomes proposals from external investigators interested in 
collaboratively using the dataset while it is still proprietary, by 
contacting the team leaders at gama@gama-survey . org. 
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